Topological defects in graphene: dislocations and grain boundaries 
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Topological defects in graphene, dislocations and grain boundaries, are still not well understood 
despites the considerable number of experimental observations. We introduce a general approach 
for constructing dislocations in graphene characterized by arbitrary Burgers vectors as well as grain 
boundaries, covering the whole range of possible misorientation angles. By using ab initio calcula- 
tions we investigate thermodynamic and electronic properties of these topological defects, finding 
energetically favorable symmetric large-angle grain boundaries, strong tendency towards out-of- 
plane deformation in the small-angle regimes, and pronounced effects on the electronic structure. 
The present results show that dislocations and grain boundaries are important intrinsic defects in 
graphene which may be used for engineering graphene-based nanomaterials and functional devices. 

PACS numbers: 61.48.Gh, 73.22.Pr, 61.72.Lk, 61. 72. Mm 



I. INTRODUCTION 

The isolation of graphene, a two-dimensional (2D) 
material with extraordinary physical properties, has 
opened new horizons for physics exploration and future 
technology In 2D, properties of materials can be heav- 
ily affected by structural irregularities. Graphene edges 
and point defects such as vacancies have been extensively 
investigated over the past few years. However, these 
types of disorder have to be distinguished from dislo- 
cations and grain boundaries, structural defects charac- 
terized by the finite values of their respective topological 
invariants. Burgers vectors and misorientation anglesi^ 
Such topological defects introduce non-local disorder into 
the crystalline lattice. Surprisingly, dislocations and 
grain boundaries in graphene are still not well under- 
stood despite the growing number of experimental obser- 
vations. 

The first experimental results date back to the scan- 
ning tunneling microscopy (STM) studies of tilt grain 
boundaries on graphite surfaces,^ fueled by their con- 
fusion with biological macromolecules<^ii More recently 
it has been shown that grain boundary defects have a 
dramatic influence on the local electronic properties of 
graphite. ^'^ An individual dislocation in free-standing 
graphene layers has been imaged using transmission elec- 
tron microscopy (TEM).^^ Topological defects resulting 
from either kinetic factors or substrate imperfections 
have also been reported for epitaxial graphene grown on 
SiC,— Ir(lll)i2iH and polycrystalHne Ni surfaces-^^ 

Here, we describe a systematic approach for con- 
structing arbitrary dislocations and grain boundaries in 
graphene starting from disclinations as the elementary 
topological defects. Then, by using ab initio calcula- 
tions we explore energetic and electronic properties of 
the proposed structures finding a number of intriguing 
features such as two energetically favorable symmetric 
large-angle grain boundaries, strong tendency towards 
out-of-plane deformation in the small-angle regimes, and 
pronounced effects on the electronic structure. Our re- 



sults highlight the possible important role of dislocations 
and grain boundaries in practical graphene samples. 

The present paper is organized in the following man- 
ner. In Section [TT] we describe our first- principles compu- 
tational methodology. Section llll Al presents a systematic 
approach for constructing atomic structures of disloca- 
tions and grain boundaries in graphene. Sections IIII Bl 
and IIII CI are devoted to the discussion of energetics and 
electronic structure of the constructed topological de- 
fects, respectively. Section llVl concludes our work. 



II. COMPUTATIONAL METHODS 

First-principles calculations have been performed us- 
ing the spin-polarized density functional theory (DFT) 
scheme implemented in the SIESTA codei^^ The general- 
ized gradient approximation (GGA) exchange-correlation 
density functional^— was employed together with a 
double-(^ plus polarization basis set, norm-conserving 
pseudopotentials^'' and a mesh cutoff of 200 Ry. 
The computational model involved two parallel equally 
spaced grain boundaries in a rectangular simulation su- 
percell in order to satisfy periodic boundary conditions 
(see Fig. [T|) . The distance between the neighboring dis- 
locations along the boundary line, and thus the misori- 
entation angle 6, are changed by varying the dy supercell 
dimension. The supercell dimension was ~4 nm in all 
studied models. The chosen supercell construction allows 
one to reduce the error due to elastic interactions between 
the neighboring grain boundariesJ^ We verified that a 
larger inter-boundary separation (dx = 8 nm) produces 
only a negligible change in the calculated grain boundary 
energy in both small-angle and large-angle regimes. Both 
atomic coordinates and supercell dimensions were opti- 
mized using the conjugate-gradient algorithm and a 0.04 
eV/A maximum force convergence criterion. The Bril- 
louin zone was sampled using 2 fc-points along the x axis 
and a consistent number of approximately 8/dy fc-points 
along the y axis {dy in nm). The scanning tunneling 
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FIG. 1: (a) Schematic illustration of the dx x dy rectangu- 
lar simulation supercell with two dislocations (filled symbols) 
separated hy dx/2. The periodic images of the dislocations are 
shown as empty symbols. The dashed lines depict the grain 
boundary lines, (b) One of grain boundary models {0 = 9.4°) 
used in the present study. The simulation supercell is indi- 
cated. 



microscopy (STM) images were simulated using a pre- 
viously developed metho d^^i^° employing the calculated 
local density-of-states in the energy window of ±0.6 eV 
around the charge neutrality point. 



III. DISCUSSION OF RESULTS 

A. Atomic structure of topological defects: 
a systematic approach 

In truly 2D materials only edge dislocations are pos- 
sible since the Burgers vector &, a topological invariant 
which reflects the magnitude and direction of the crys- 
talline lattice distortion produced by a dislocation, is con- 
strained to lie in the material's plane. One can imagine 
such dislocation as a result of embedding a semi-infinite 
strip of width |6| into an otherwise perfect 2D crystalline 
lattice. As a guiding rule for constructing atomic struc- 
tures of dislocations in graphene we assume that the dis- 
location core is free from under- or over-coordinated car- 
bon atoms; that is, we aim at minimizing the energy of 
the dislocation core and, thus, the total formation en- 
ergy of the dislocation.— To develop such construction 
we adopt a membrane theory approach which views a dis- 
location as a pair of positive and negative disclinations, 
i.e. topological defects obtained by removing and adding 
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FIG. 2: (a) Positive (s = 60°) and negative (s = -60°) 
disclinations in graphene are produced by either removing or 
adding a 60° wedge (shaded area) of material without chang- 
ing the coordination of carbon atoms, (b) A pair of comple- 
mentary disclinations is equivalent to a dislocation: a negative 
disclination inserts a 60° wedge while a positive disclination 
removes such a wedge within one of the seven equivalent sec- 
tors. The introduced amount of material (shaded area) can 
also be viewed as a semi- infinite strip of width 



a semi-infinite wedge of material to an otherwise per- 
fect crystalline lattice, respectively.^^ As shown in Fig- 
ure [^a), s — 60° (s — —60°) disclination in graphene 
contains a five (seven) membered ring in its core while 
the original three-fold coordination of all carbon atoms 
is preserved. Figure [^b) schematically shows the equiv- 
alence of a pair of complementary disclinations to a dis- 
location. Moreover, we find that on graphene lattice the 
distance between two disclinations, |a|, is related to the 
resulting Burgers vector 6 by a simple relation, |a| = \b\ 
(for proof see Appendix |^ . 

Since any Burgers vector 6 is a proper translational 
vector of graphene lattice, i.e. b = ndi + ma2 (ai,2 = 
(3(icc/2, ±A/3dcc/2); dec = 1-42 A, the nearest neigh- 
bor interatomic distance in graphene), we will use the 
pair of integers (n, m) as a descriptor of dislocations in 
graphene. This notation is analogous to the chirality 
indices used to describe the structure of carbon nan- 
otubes. The core of the shortest Burgers vector dis- 
location (1,0) (|6(i,o)l = VSdcc = 2.46 A) contains an 
edge-sharing heptagon-pentagon pair as shown in Fig- 
ure [3fa). The (1,0) dislocation inserts a semi-infinite 
strip of atoms along the armchair high-symmetry di- 
rection in graphene while its Burgers vector is oriented 
along the zigzag direction. This simplest dislocation 
structure has been extensively studied in the context 
of plastic deformation of carbon nanotubes^ nanotubes 
junctions^ as well as graphene itselfi ^^i^^ The second 
member of the family, the (1,1) dislocation, has a larger 
Burgers vector (|6| = Sdcc — 4.23 A) and inserts a semi- 
infinite strip along the zigzag direction of graphene (see 
Fig-Efb)). Alternatively, the core of the dislocation with 
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FIG. 3: (a-c) Atomic structures of (1,0) and (1,1) dislocations, and a (1,0) + (0,1) dislocation pair, respectively. The dashed lines 
delimit the introduced semi-infinite strips of graphene originating at the dislocation core. Non-6-membered rings are shaded. 
(d,e) Atomic structures of the 6 — 21.8° (LAGB I) and the 6 = 32.2° (LAGB II) symmetric large-angle grain boundaries, 
respectively. The dashed lines show the boundary lines and the solid lines definite angles 6i and 62- (f) Buckling of the 
graphene layer due the presence of a (1,0) dislocation. 



the same Burgers vector can be constructed from two 
1^(1,0)1 = 2.46 A dislocations, (1,0) and (0,1), e.g. as 
shown in Figure [2{c). The simple method outlined above 
can be used to build dislocation with even longer Burgers 
vectors, inevitably leading to larger elastic energies. 

Grain boundaries, the interfaces between the domains 
of material with different crystallographic orientations, 
are commonly viewed as periodic arrays of dislocations.^^ 
Particularly, in 2D materials such as graphene, one- 
dimensional (ID) chains of edge dislocations constitute 
tilt grain boundaries. Mutual orientation of the two crys- 
talline domains is described by the misorientation angle 
e ^ 01+ 02 {9 e (0°,60°) in graphene), a topological 
invariant defined as shown in Fig. ^d). Another pa- 
rameter tp — \di — E (O°,0) describes the inclina- 
tion of the boundary line with respect to the symmet- 
ric configuration = 0°). We limit our consideration 
to only symmetric ones since asymmetric configurations 
tend to result in diverging elastic energiesj^^ Importantly, 
due to the presence of two high-symmetry directions in 
graphene, armchair and zigzag, both misorientation an- 
gles close to 0° and 60° can be considered as small-angle 
grain boundaries along these two directions, respectively. 
Aligning (1,0) dislocations along the grain boundary line 
results in a discrete set of misorientation angles in ac- 



cordance with Frank's equation^! 

= 2arcsinl^, (1) 
2"(i,o) 

where d(i,o) is one of the possible values for the distance 
between the neighboring dislocations. Large values of 
0) correspond to small-angle grain boundaries along 
the armchair direction. The closest possible packing of 
(1,0) dislocations results in the large-angle grain bound- 
ary structure shown in Fig. ^d). This configuration char- 
acterized by — 21.8° (LAGB I) has already been sug- 
gested in the literature In order to cover the the range 
of between 21.8° and 60°, it is necessary to introduce 
another type of dislocations, e.g. (1,1) dislocations: 

6* = 60° - 2 arcsin l^^i^ . (2) 
2"(i,i) 

The smallest value of gives rise to the LAGB I 

structure rotated by 180°. Large separations 
correspond to small-angle grain boundaries along the 
zigzag direction. Alternatively, small-angle grain bound- 
aries along this direction can be constructed using the 
(1,0)-|-(0,1) pairs with the densest possible packing of dis- 
locations leading to the structure LAGB II with — 32.2° 
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(Fig. [3je)). Hence, it is possible to construct symmetric 
grain boundaries covering the whole range of 9 from 0° 
to 60° by using (1,0) dislocations and either (1,1) dislo- 
cations or (1,0) + (0,1) dislocation pairs. We stress that 
the present construction is equally applicable to tilt grain 
boundaries in both graphene and graphite. 



B. Energetics of topological defects 

In order to determine the energetically preferred struc- 
tures and to understand the basic thermodynamic prop- 
erties of grain boundaries in graphene, we perform first- 
principles calculations on models containing a pair of 
complementary dislocations in periodic 2D supercell as 
described in Section [TTl The results are presented as a 
diagram of grain boundary energies per unit length 7 
as a function of 9 (Fig. H]). We first discuss the case 
of perfectly flat grain boundaries (filled symbols) which 
corresponds to the limit of strong binding to a fiat sub- 
strate and analogous to the case of grain boundaries in 
bulk materials. The diagram clearly reveals both arm- 
chair and zigzag small-angle regimes with grain bound- 
ary formation energies converging to zero for — >■ 0° and 
— >■ 60°. A detailed study of the armchair small-angle 
region [0 < 10°) shows that the grain boundary energies 
are well described by the Read-Shockley equation^i 



1(0') - 



47r(l - 1^) 



9'(l + ln 



\b\ 
27rro 



ln9'), (3) 



where fj, is the shear modulus, i' the Poisson's ratio. The 
core radius rg encompasses the energy of the disloca- 
tion core. In Eq. ©, 0' = 6* or 9' = 60° - 9 for arm- 
chair and zigzag small-angle grain boundaries, respec- 
tively. Using the values of elastic constants which cor- 
respond to our first principles model of graphene (see 
Appendix [B| a least-squares fit to the Read-Shockley 
equation (solid red curve in Fig. 2]) yields tq = 1.2 A. 
This value is in good agreement with the recently re- 
ported ro — 0.96 A fitted to local density approximation 
calculations As we outlined above, there are several 
grain boundary structures possible for 6* > 21.8°. In or- 
der to determine the lowest energy structure, we compare 
the energies of grain boundaries constructed from (1,1) 
dislocations and (1,0) + (0,1) dislocation pairs. In addi- 
tion, for 9 > 42.1° the grain boundary can be constructed 
either from equally spaced (1,0) and (0,1) dislocations 
(disperse case) or from closely bound pairs (paired case). 
Figure|l]shows that the disperse (1,0) + (0,1) grain bound- 
aries in flat graphene are the lowest-energy structures for 
9 > 42.1°. More generally, this also implies that only 
the shortest Burgers vector (1,0) dislocation is sufficient 
for constructing the most stable grain boundary struc- 
tures at any given 9. Remarkably, the two large-angle 
structures discussed above, LAGB I and LAGB II, have 
particularly low formation energies of 0.338 eV/A and 
0.284 eV/A, respectively. Favorable energetics suggests 
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FIG. 4: (color online). Grain boundary energy per unit 
length 7 as a function of misorientation angle 9 for various flat 
(filled symbols) and buckled (open symbols) grain boundary 
structures. The two energetically favorable large-angle grain 
boundaries, LAGB I and LAGB II, are labeled. Solid curve 
shows the Read-Shockley equation fit (ro = 1.2 A) for the flat 
small-angle armchair grain boundaries. Dashed curve shows 
the asymptotic linear dependence of 7 for the buckled small- 
angle armchair grain boundaries {Ef — 7.5 eV). 



possible abundance of these two structural motifs. More- 
over, for all possible values of 9 the grain boundary en- 
ergies are well below the energies of ~1 eV/A predicted 
for graphene edges <^ 

The case of free-standing 2D materials is notably dif- 
ferent since buckling in the third dimension allows an ex- 
change of in-plane elastic energy for bending energy. This 
leads to efficient screening of the in-plane strain field re- 
sulting in the finite formation energies of dislocations^— 
While we find that the large-angle grain boundaries in 
graphene are flat, for 9 < 21.8° and 9 > 38.2° buckling 
effectively reduces the grain boundary energies (Fig. |4l 
empty symbols). In the small-angle regimes, grain 
boundary energy is expected to scale linearly with 9': 



(4) 



where Ef is the formation energy of the dislocation. By 
extrapolating the grain boundary energies to = 0°, 
we obtain a formation energy of 7.5 eV for the (1,0) 
dislocation (see Appendix [C]). This value is compara- 
ble to the formation energies of typical point defects in 
graphene, e.g. vacancies (7.6 eV) and Stone- Wales de- 
fects (4.8 eV).'^^ Out-of-plane buckling results in a prolate 
hillock appearance of dislocations on the flat graphene 
surface (Fig. [Sljf)) in agreement with the experimental 
observations of Coraux et al}^ The height of the protru- 
sion around the (1,0) dislocation is ~3 A and its top is 
shifted with respect to the dislocation core. Interestingly, 
out-of-plane distortion makes the paired case (1,0)-|-(0,1) 
grain boundaries more stable, thus inverting the sign of 
effective interaction between the dislocation dipoles. The 
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FIG. 5: (color online), (a) Calculated density-of-states plots 
for the small-angle armchair {9 = 7.3°) and the possible con- 
figurations of small-angle zigzag {6 = 49.5°) grain bound- 
aries. The plots correspond to the values averaged over 2 nm 
wide interface regions. The dotted line shows the density-of- 
states of the ideal graphene. (b) Calculated density-of-states 
plots for the large-angle grain boundary structures LAGB I 
and LAGB II. (c) Simulated STM images of the individual 
dislocations in small-angle grain boundaries and large-angle 
structures. The images cover 2 nm x 2 nm areas. 




FIG. 6: The calculated band structures for the models of 
(a) LAGB I and (b) LAGB II large-angle grain boundaries 
along at fcx = 0. Correspondence between the van Hove 
singularities in Figure [5jb) and the band extrema in the band 
structure plots is highlighted with arrows. 



further corroborated by considering the electronic band 
structures of large-angle grain boundary models shown 
in Figure |B1 However, we do not observe any zero- 
energy states or defect-induced magnetic moments typ- 
ical of zigzag edges^i^i and single-atom defects^"— in 
graphene. Pronounced changes in the low-energy part of 
electron spectrum make it possible to identify the dis- 
cussed extended defects using STM. In order to facilitate 
the attribution of experimental observations to the pro- 
posed structures we provide their simulated atomic-scale 
STM fingerprints (Fig. [SJc)). The common feature of 
all images is the crescent- or ring-shaped appearance of 
5-membered rings. The STM images of grain bound- 
aries formed by (l,0)-f (0,1) dislocations show the lack of 
mirror symmetry compared to the (1,0) or (1,1) derived 
structures, as follows from the atomic structures of these 
grain boundaries. 



situations in which graphene is bound to substrate can 
be viewed as intermediate between the flat and buckled 
regimes. 



C. Electronic structure of topological defects 

Finally, we address the electronic structure of topolog- 
ical defects in graphene. Figures [SJa) and EJb) show the 
calculated density-of-states plots for small-angle and sta- 
ble large-angle grain boundaries in graphene. All stud- 
ied defect configurations introduce van Hove singular- 
ities within 0.5 eV below and above the Dirac point 
(Ep = eV), in accordance with the scanning tun- 
neling spectroscopy (STS) observations for the major- 
ity of grain boundaries in graphite reported in Refs. 
The van Hove singularities are the signatures of one- 
dimensional states localized at the interface as shown 
by simulated STM images (Fig. EJc)). This relation is 



IV. CONCLUSIONS 

We have developed a systematic approach for con- 
structing atomic structures of topological defects in 
graphene. Our first-principles calculations revealed a 
number of intriguing features in the energetics of grain 
boundaries. In particular, we have found two large-angle 
grain boundary structures with particularly low forma- 
tion energies as well as two distinct small-angle regimes 
which correspond to the grain boundaries oriented close 
to the armchair and zigzag directions, respectively. In 
free-standing graphene the small-angle grain boundaries 
show pronounced tendency to an out-of-plane buckling 
which further reduces their formation energies. We have 
also found that all the studied topological defects have 
strong effects on the electronic structure and can be iden- 
tified using STM. These results show that dislocations 
and grain boundaries are important intrinsic defects in 
graphene which may be used for engineering graphene- 
based nanomaterials and functional devices. 
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Appendix A: Relation between the dislocation 
dipole a and the Burgers vector b 

Theorem. // a dislocation is constructed from a pair 
of s = ±7r/3 disclinations, then its Burgers vector b and 
vector a connecting the disclinations are related as 



\a\ = \b\. 



(Al) 



Proof. We consider the following construction in which 
a negative s = — 7r/3 disclination inserts sector AHB at 
point H and a complementary positive s — t: /i disclina- 
tion removes sector A" PB" at point P in a continuous 
two-dimensional sheet. 

The condition of continuity requires AH = BH 
AB), WH = WH as well as WW = WH. In addition, 
we require that, after the described procedure, pairs of 
segments A' A" and B'B" as well as AP and BP form 
straight lines. That is. 



Hence, 



A' A" + B'B" =AP + BP^ \b\. 



APA" + BPB" = TT. 



(A2) 



(A3) 




By construction AHB = A"PB" = 7r/3, thus APB = 

27r/3 a.nd'APB + AHB = IIAP + HBP = tt. Due to the 
latter property, quadrilateral AH BP can be inscribed in 
a circle and, thus, from Ptolemy's theorem it follows that 



Finally, 



AB ■ HP = BH ■ AP + AH ■ BP. (A4) 



\a\^ HP = AP + BP = \b\. (A5) 



Appendix B: Structural and elastic constants of 
graphene from first principles 

In order to fit the results of our calculations of small- 
angle grain boundaries to the continuum-model Rcad- 
Shockley equation, we use elastic constants and the in- 
teratomic distance of graphene which correspond to the 
present first-principles model of graphene. Elastic con- 
stants are obtained from the constrained variable cell cal- 
culations in which one of the rectangular supercell di- 
mensions was fixed while the other one varies in order to 
minimize the total energy. This allowed us to determine 
Young's modulus E and Poisson's ratio v while the shear 
modulus /i was calculated using the following relation)^ 



E 



2(1 + 



(Bl) 



The calculated moduli correspond to 3.35 A thickness of 
the graphene layer. Our values are in good agreement 
with other values reported in literature (see Table H]). 



Appendix C: Formation energy of the buckled (1,0) 
dislocation in graphene 

Figure S] shows that the energies of the buckled small- 
angle armchair grain boundaries do not achieve the ex- 
pected linear dependence (Eq. ([4])) in the range of stud- 
ied misorientation angles 9. The maximum separation 
D between the neighboring (1,0) dislocations along the 
boundary line we could afford in our demanding first- 
principles calculations is ^4 nm (corresponds to — 



A)4— 1 — \ 

FIG. 7: Sketch of the geometric construction of a dislocation 
in a continuous sheet. Negative and positive 7r/3 disclinations 
are placed at points H and P, respectively. 



TABLE I: The values of interatomic distance dec, Young's 
modulus E, Poisson's ratio u and the shear modulus /j, of 
graphene obtained from the present first-principles calcula- 
tions. The results are compared to values reported in litera- 
ture. 

dec (A) E (TPa) u n (GPa) 

this work 033 052 0206 436 

Exp. (graphite)^ 1.02±0.03 0.16±0.03 

Exp. (graphene)^ l.OiO.l 
Theory^J. 1.050 0.186 

Theory^? 1.01±0.03 0.21±0.01 
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3.5°). At this D, the screening of the in-plane elastic 
field is still insufficient to decouple the neighboring dis- 
locations along the grain boundary direction. However, 
we observe that the energy per dislocation jD shows a 
clear linear dependence with 1/D. By extrapolating the 
vahics of grain boimdary energies for < 10° to the limit 
of 1/D = (that is, = 0°), we obtain an estimate of the 
formation energy of an isolated buckled (1,0) dislocation 
Ef = 7.5 eV. 



FIG. 8: Grain boundary energies per dislocation 'yD as a 
function of inverse distance D between the neighboring dislo- 
cations. Dashed line shows the least-squares fit. 
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